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We study non-adiabatic corrections to multibaryon systems within the bound state approach to 
the SU(3) Skyrme model. We use approximate ansatze for the static background fields based on 
rational maps which have the same symmetries of the exact solutions. To determine the explicit 
form of the collective Hamiltonians and wave functions we only make use of these symmetries. Thus, 
the expressions obtained are also valid in the exact case. On the other hand, the inertia parameters 
and hyperfine splitting constants we calculate do depend on the detailed form of the ansatze and 
are, therefore, approximate. Using these values we compute the low lying spectra of multibaryons 
with B < 9 and strangeness 0,-1 and —B. Finally, we show that the non-adiabatic corrections do 
not affect the stability of the tetralambda and heptalambda found in a previous work. 
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I. INTRODUCTION 
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In the last few years there have been several important developments in the determination of the lowest energy 

^■—2 ■ skyrmion configurations pHl . This type of solutions are essential for the understanding of multibaryons and, perhaps, 

nuclei in the framework of the topological chiral soliton models. So far, these models have proven to be useful for the 

description of quantities such as the masses, strong and electro-magnetic properties of the octet and decuplet baryons, 

■?G baryon-baryon interactions, etc. (see e.g. Refs. plpf and references therein). The knowledge of the properties of 

£~H. the multiskyrmion configurations opens the possibility of studying more complex baryonic objects. In fact, several 

O ,' investigations concerning non-strange multiskyrmion systems have been reported in the literature (see, e.g., Refs. 

[|6| |l0|). Of particular interest are, however, the strange multibaryons. Perhaps the most celebrated example is the 

*G H dibaryon predicted in the context of the MIT bag model more than twenty years ago pT[ . This exotic has been 

studied in various other models, including the Skyrme model fl2}-[l5||, but its existence remains controversial both 

theoretically and experimentally. It has also been speculated that strange matter could be stable jjl|. This has lead 

to numerous investigations of the properties of strange matter in bulk and in finite lumps (for a recent review see Ref. 



Jr7|). Moreover, with the new heavy ion colliders there is now the possibility of producing strange multibaryons in 
the laboratory |l8|]. In this situation the study of multibaryon systems within the SU(3) Skyrme model appears to be 
quite interesting. A first step in this direction has been reported in Ref. [fill where the rational map approximation 
pO | to the multiskymion fields was used to describe the multibaryon configurations within the bound state approach 
pi | to the ST/ (3) Skyrme model. Within this approach strange (multi)baryons appear as systems of kaons bound to 
a background skyrmion configuration. To find the kaon binding energy one has to solve the corresponding eigenvalue 
problem. For a general background this is a very hard numerical task since one has to deal with several couple partial 
differential equations. However, this problem is greatly simplified if one introduces the (approximate) rational maps 
ansatze for the multiskyrmion configurations. The construction of these ansatze is based on the analogy between 
monopoles and skyrmions and requires that the approximate solutions have the same symmetries than the exact 
numerical solutions. In fact, it is now known that up to B = 9 these configurations are very symmetric. Namely, 
for B = 2 the solution corresponds to an axially symmetry torus while configurations with B = 3 — 9 possess the 
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symmetries of the platonic polyhedra. In contrast with the exact solution, however, the rational map approximation 
assumes that the modulus of the static pionic field is radially symmetric while its direction depends only on the polar 
coordinates. It was shown in Ref. |2(J that this represents a very good approximation. Once the rational maps are 
introduced the kaon eigenvalue problem reduces, for each baryon number, to one radial eigenvalue equation. The 
corresponding results have been given in Ref. p9[ . In such reference, however, non-adiabatic effects were neglected. 
These effects appear when one performs the collective quantization of the system. It should be stressed that it is only 
at this stage when the spin and isospin quantum numbers are well defined and splitting between the corresponding 
states appears. The purpose of the present work is to carry out the collective quantization of the bound multisoliton- 
kaon systems. This requires to pay special attention to their symmetries which impose severe constraints on the 
possible quantum numbers and wave functions. 

This paper is organized as follows. In Sec. II we provide a brief description of the model with special emphasis on 
the effect of the non-adiabatic corrections. In Sec. Ill we describe in detail how to obtain the collective Hamiltonian 
for the different baryon numbers, while in Sec. IV we focus on the corresponding wavefunctions. It should be 
noticed that since the discussions in these two sections rely only on the symmetries of multiskyrmion configuration 
the corresponding results hold true also for the exact solutions. In Sec. V we present the numerical results and in 
Sec. VI our conclusions. Finally, in the Appendix we give the explicit form of the rational maps used in the present 
work. 



II. THE MODEL 



We start with the effective action of the SU(3) Skyrme model supplemented with an appropriate symmetry breaking 
term ||. Expressed in terms of the 5?/(3)-valued chiral field U{x) it reads 



r= fd 4 x ^Tr[d^Ud^] + ^Tr^d^U^d.U} 2 ]] +T WZ +Tsb 



(1) 



where / w is the pion decay constant ( = 93 MeV empirically) and e is the so-called Skyrme parameter. In Eq. ([[]), 
the symmetry breaking term T$b accounts for the different masses and decay constants of the pion and kaon fields 
while Twz is the usual Wess-Zumino action. Their explicit forms are 
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where A 8 is the eighth Gell-Mann matrix and m^ and thk represent the pion and kaon masses, respectively, and fx 
is the kaon decay constant. 

We proceed by introducing the Callan-Klebanov ansatz for the chiral field |El[ 

(4) 



In this ansatz, Uk is the field that carries the strangeness. Its form is 



Uk — exp 



.a/2 / K 
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where K is the usual kaon isodoublet K = 



K+ 



,-,: j . The other component, U ni is the soliton background field. It 
is a direct extension to 5/7(3) of the SU{2) field, i.e., 
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Replacing the ansatz Eq.(H) in the effective action Eq.(pl) and expanding up to second order in the kaon fields we 
obtain the Lagrangian density for the kaon-soliton system. In the spirit of the bound state approach this coupled 



system is solved by finding first the soliton background configuration. For this purpose we introduce the rational map 
ansatze 120(1 



with 
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where we have assumed that F = F(r), and R = R(z) is the rational map corresponding to winding number B = n. 
Here, r is the usual spherical radial coordinate whereas the complex variable z is related to the other two spherical 
coordinates (0, 4>) via stereographic projection, namely, z = tan(0/2) exp(z^). The resulting expression for the soliton 
mass per unit baryon is (in what follows s = sin F; c = cos F) 
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The profile function F(r) is obtained by minimizing M so i subject to the boundary conditions F(0) = 7r and F(oo) = 0. 
In using these boundary conditions we are assuming that all the extra winding number is obtained from the angular 
dependence of w. The angular integral T is 



r r " f^nm- ^m 2 l fjidzdz_ f l + \z\ 2 
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In order to find the lowest soliton- kaon bound state we write the kaon field as fi4J , [l5| , 

K rAr,t) = k(r,t) T-n XT; , 

where x is a two-component spinor. 

The diagonalization of the corresponding kaon Hamiltonian leads to the eigenvalue equation 



-\d r (r 2 h n d r ) +m 2 K + V* ff - /„e 2 - 2 A„ e„ 
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Details on how to obtain this equation as well as the explicit expression of the radial functions /„ , h n , X n and V n can 
be found in Ref. [§§. 

To obtain the hyperfine corrections to the multibaryons masses we proceed with the semi-classical collective coordi- 
nates quantization method, where the isospin and spatial rotations are treated as the zero modes. Then, we introduce 
the time-dependent spatial rotations R and the isospin rotations A such that 
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R Au„ A~ 



RAK. 

The angular velocities with respect to the body fixed frame are given by 

\R~ R) = tabc^c , 
\ J ab 

A^A^-f-uJ. 
2 

Replacing in the effective action we get the collective Lagrangian 

Lcou = -m so1 + - [e( b n a n b + ei b u a u b + 2 e^ fi a ^ 6 ] - (c( b n a + c4w a ) r b , 

where a,b — 1,2,3 and T b = (x^ T bX)/2 is the kaon spin. 

The moments of inertia <d ab and hyperfine splitting constants c ab appearing in Eq. (|l7|) are given by 
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where the radial integrals mi, m.2, rfi and c?2 are 
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and the angular integrals 
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The numerical values of these angular integrals depend only of the particular form of the ansatz for h and not on the 
detailed form of the effective action and its parameters. For the rational maps listed in the Appendix all the matrices 
Eqs. (p7|)-(|32|) are diagonal. As we shall see in the next section, this is a direct consequence of the symmetries of these 
ansatze. The corresponding values of the diagonal elements are listed in Table I. Note that when all the diagonal 
elements are equal we list just one. Also listed in Table I are the values of X. 

Given L co u, the canonical momenta are then defined in the usual way 
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where we have used that, for the cases we are interested in, all the inertia and hyperfine splitting constants are 
diagonal and thus denoted with a subindex a — 1,2,3 the corresponding diagonal elements. Depending on whether 
A a = 0^0^ — ( < d a vI ) 2 vanishes or not we have to follow a somewhat different procedure to obtain the collective 
Hamiltonian. We consider first the case in which A a ^ for all values of a. In this case the relations Eq.(33) can be 
inverted and the collective Hamiltonian results 
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where 



and 
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If there exist, however, some values i for which Aj = there appears a relation between J^, J^ and T^. It reads 
Using this relation it is not difficult to show that the collective Hamiltonian becomes 



jj CO U = ^ jjaaU + J- ^/" . (39) 
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and the total multibaryon mass results 

M = n M soi + \S\e n + E rot (40) 

where S is the multibaryon strangeness and E rot the expectation value of H rot in the corresponding wavcfunction. 
In the next section we will determine the precise form of the collective Hamiltonians for each baryon number. 

III. COLLECTIVE HAMILTONIANS 

The minimum energy multiskyrmion configurations are symmetric under certain groups of transformations . With 
the exception of the B = 1 and B = 2 cases where these symmetry groups are continuous (0(3) and -Doo/i, respectively), 
these transformation groups have a finite number of elements. In this section we will see how the symmetries of the 
multiskyrmion configurations impose severe constraints on the detailed form of the collective Hamiltonian. For the 
B < 4 cases this has already been discussed in the literature using various arguments. Here, we will extend such 
analysis within a unified framework. It is important to notice that all the discussions and results that follow are based 
only on the symmetries of the multiskyrmions. Therefore, they will hold not only for the approximate configurations 
based on the rational maps but also for the exact ones obtained from numerical minimization. 

The task here is to determine the precise structure of the inertia and hyperfine splitting tensors, namely, which 
elements of those tensors vanish and how many of the remaining non-zero elements are independent for each baryon 
number. First, we note that each operation of the abstract group G is represented by a pair of operations {g, D g } which 
act in spin and isospin spaces, respectively. The pion field in Eq.(ph is invariant under these combined operations, 

f- 7 ?(f)=I? s T- 7 r( 5 - 1 f)(^ 9 ) t . (41) 

Given the form used for the kaon fieldR, Eq.(R), this invariance implies that the action of the group element on the 
kaon field is also represented by D g . In fact, 

D g K Tz (r,t) =K Tz (gr,t), (42) 

which means that the symmetry operation acting on the kaon field is just given by the representation of the isospin 
operation D g in the T-space. Thus, the I and T operators transform in the same way under elements of G. This 
shows that it is enough to perform the explicit analysis only for the inertia tensors. Once this is done the results for 
the hyperfine splitting constants can be easily obtained noting that in Eq.(pfl), c J ab plays a role similar to that of K^[, 
while c 7 fc to that of K^ b . 



x This ansatz can be easily generalized if the exact numerical soliton configuration is used instead of the approximation based 
on rational maps. 



The inertia tensors can be diagonalized by an appropriate choice of the spatial and internal reference frames, 
and this is in fact what happens for the rational map ansatze given in the Appendix. Consider first the case for 
the spin. The spin generators J a transform under G in some (possibly reducible) representation. The number of 
independent diagonal components of the inertia tensor (moments of inertia) will be equal to the number of irreducible 
representationsn (irreps) of G into which this representation breaks, since the combination K^ b J a Ji must be a scalar 
under G. The spin generators belong to the 1 + irrep of 0(3) which for the cases we will consider below breaks 
into either a 3-dim irrep or as the sum of 1- and 2-dim irreps of G. In the first case there is only one moment of 
inertia and the spin Hamiltonian is proportional to ^2 a J a J a , while in the second case there are two moments, and 
the Hamiltonian contains the terms Jf + 3\ and J 2 . The same argument holds for the other collective operators. 

An important remark is the following. While there is a one-to-one correspondence between g and the elements of G, 
this is not necessarily the case for the operations D g . In other words, it could happen that the same D g is associated 
with two (or more) different elements in spin space. In this case, the operations D g do not span the full group G 
but a subgroup of it. As a consequence, the generators J a and I a (T a ) could transform in different representations of 
G. This would imply that the corresponding mixing inertia would vanish. Below we see that this happens for some 
values of B. 

Let us consider now the multiskyrmion configurations case by case. The B — 1 skyrmion is spherically symmetric 
g. Thus, the relevant symmetry group G is 0(3). In this case, g = D g and both J and / are in the 3-dim irrep 1 + . 
Using the arguments given above we have 

ei = e J , e£ = e', e? = e M , c J a =c J , <£ = <?. (43) 

Since in this case we are dealing with a continuous group the equality between the representation of the group elements 
in spin and isospin spaces can be written in terms of corresponding generators of the algebra. Namely, we obtain the 
relation J a = I a + T a . From Eq.(|33]) this implies 

e J = e 7 = e M , c 1 = 1 - C J , (44) 

which leads to A Q = for all values of a. Then, the collective Hamiltonian takes the well-known form 

— ( I 2 + c 2 T 2 + 2 c f ■ I\ . (45) 
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As already mentioned, the B = 2 lowest energy skyrmion configuration has the symmetry of a torus |l[ which 
implies G = D^h. Choosing the symmetry axis along the z-direction we obtain that the third components of the 
momenta are in the 1-dim S~ while the other two components arc in the 2-dim irrep n g . Since rotations along the 
z-axis form a continuous subgroup of -Dqo/i we obtain for the terms containing third components of the momenta a 
result similar to that of B = 1 , 

e 3 J = ®i = ef , 4 = 1 - 4 , (46) 

which leads to A3 = 0. For the other components Ai^ 7^ since the Ci along those axes only form finite subgroups 
of G. Consequently, the corresponding component of the different type of inertia and splitting constants need not to 
be equal and the B = 2 collective Hamiltonian reads 

H C B % = K{ (J 2 - J 2 ) + K[ {I 2 It) + KM)* (T 2 - T 2 ) 



+K[ci(I + T_ + I_T + )+ i -^±^l 
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For the rest of the baryon numbers under consideration, B = 3 — 9, the symmetry group G is finite P,p). Therefore, 
A a never vanishes for all those baryon numbers and the collective Hamiltonian will have the general form Eq.(p5[). 
There can be, however, some further simplifications depending on the way in which the symmetry is realized in spin 
and isospin spaces. 

The symmetry group of the B — 3 solution is G = Td- In this case, we have that g = D g for all the elements of G 
[|7J. Thus, the components J a , I a and T a are in the 3-dim irrep F%. The collective Hamiltonian reads 



2 The character tables containing the list of irreps of the groups we are interested in can be found, e.g., in Refs. |22| and 
We follow the conventions of Ref. |22| . 



H B % = K J J 2 + K 1 I 2 - 2K M I • J + 2K J c J J- f + 2K I c I I ■ T 

+ K* &-&")* <*' (C " J)2 + K ' (C_/)2 + 2RM Cl ^ ^ ' (48) 

In the case of B = 4 the relevant symmetry group is Oh- As discussed in Ref. pf, for the minimum energy 
configuration this symmetry is realized in such a way that the elements D g cover four times the D 3 d subgroup. As 
a result, I\ (Ti) and Ii (T2) are in the 2-dim irrep E gi I 3 (T 3 ) in the A 2g irrep and the components of J lie in the 
3-dim irrep T\ g . We see then that the mixing inertia and spin splitting tensors vanish. The resulting form of the 
corresponding collective Hamiltonian is 

#s=4 = K J J 2 + K[ (1 + c[ff + {Kl - K{) I 2 

+2(Kl4 K[c[) h T 3 + (Ki(4f - K[{c[) 2 ) T 2 . (49) 

The lowest energy multiskyrmion with B = 5 has Did symmetry. In this case, there is a one-to-one correspondence 
between the realization of the group in spin and isospin spaces. It is easy to check that the third components of the 
momenta are in the A 2 irrep while the other two components in the 2-dim one E. The resulting collective Hamiltonian 

is 
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As found in Ref. ||, for B = 6 the symmetry group is D^d- Due to the way in which the generators of the group 
are realized as pairs of spin-isospin operations it is possible to show that while the spin operations cover the full D^d 
group the isospin one cover twice the Did subgroup. From the corresponding compatibility tables together with the 
compatibility table of the full rotational group we find that J3, I 3 and T 3 transform as the A2 irrep, J\ and J2 as the 
E 3 irrep and the rest as Ei irrep. Therefore, 
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The B — 7 configuration has icosahedral symmetry Ih with the symmetry realized in such a way that the components 
of the spin operators transform like the F\ g irrep while those of the isospin operators as Fi g irrep. Thus, the collective 
Hamiltonian takes the simple form 

#b=7 = K J J 2 + K 1 (I + c ! f) 2 . (52) 

For B = 8 we have to deal with the D$d group. Like the case of lower even baryon numbers the isospin operations 
do not span the full group but twice a subgroup, D 3 d in this case. We find that J 3 ,I 3 and T 3 transform as the A2 
irrep, J\ and J2 as E$ irrep and the rest as the £4 irrep. This implies that the collective Hamiltonian for B = 8 has 
the same form as the B — 6 one given in Eq.(j5lj). Finally, the B = 9 multiskyrmion has the same symmetry as the 
B = 3 one, Td- Consequently, we obtain a similar form for the corresponding collective Hamiltonian, Eq.([48|). 

IV. COLLECTIVE WAVE FUNCTIONS 

Having determined the explicit form of the collective Hamiltonian we have to find the corresponding wave functions. 
These wavefunctions have to satisfy some constraints imposed by the symmetries of the background multiskyrmion. 
For non-strange multiskyrmions this problem has been discussed by several authors p-pll. Here, we will extend such 
studies for kaon-soliton bound systems. 



The quantization of a single Skyrmion as a fermion implies that under certain symmetry operations of the classical 
multisoliton background the corresponding wave functions can pick up a nontrivial phase. These are known as 
Finkelstein-Rubinstein (FR) constraints p5| . We can generically write the constraints on the ground state as 

9 D g \g.s.) =7 9 |g-s-) , (53) 

where 7 9 = ±1 is determined according to the FR constraints. Using continuity arguments it turns out that the 
FR phases can be non-trivial only for those operations corresponding to rotations, so for our cases of interest only 
the proper subgroup of G needs to be considered. For the isospin transformations we have to take into account the 
fact that the symmetry operation also acts on the kaon field. From Eq.(|42|), however, we notice that this operation 
coincides with the one acting on the soliton isospin space. Thus, defining N = I + T, the problem basically reduces 
to that of non-strange baryons just replacing the collective isospin by N. The (proper) group generators and their 
corresponding FR phases for the configurations considered in this work were determined in Refs. ||7|Jl0f| . They are 
listed in Table || 

It is clear from Eq.(|53|) that due to the FR phases the soliton ground state might transform in a one-dimensional 
non-trivial irrep of G. Using the FR phases listed in Table O and the group character tables, the relevant 1-dim irrep 
F can be determined. We obtain that, except for the B = 5 and B = 6 cases, all the wavefunctions should transform 
as the trivial irrep of the corresponding symmetry groups. For B = 5, T is the A^ irrep of Did while for B = 6 the 
wave functions should transform as the Ai irrep of D44. 

We now need to determine the collective wavefunctions. The general procedure for arbitrary soliton backgrounds 
was discussed in p6| . We consider first the problem without strangeness. In this case we need to determine the 
functions 

I JJ Z ,II,) = J2 <*%h D i Js D U > ( 54 ) 

J3I3 

which transform under the right action of G in the irrep T of the soliton. This can be done following standard group 
theoretical methods J271]. The product representation J x I of SU{2) is in general a reducible representation of G. 
The projector operator into the irrep T is 

P r = ^E*r(5)p(.9), (55) 

11 gee 

where \G\ is the rank of the group, xr(g) the character of operation g, and p(g) the representation of g in J x / (cf. 
Eq.©) 

Pis) = DJ (9) x D^Dg) . (56) 

The eigenvalues of Pr can either vanish or be equal to one. The eigenvectors corresponding to each non-vanishing 
eigenvalue provide precisely the coefficients q^j of Eq.(J5J), and there are as many wave functions as non-zero eigen- 
values. If all eigenvalues vanish there is no collective state with the given J, /. If there is only one, the wavefunction is 
an eigenfunction of the collective Hamiltonian, and if there are more than one, the Hamiltonian has to be diagonalized 
in the subspace spanned by them. 

Let us proceed now to the case with S ^ 0. We need to find the functions!] 

\JJ Z ,II Z ,S) = J2 J / 3 l T3 Dj zj3 Dl h Kl , (57) 

which transform in irrep T under G. However, as noted above, the action of G in isospin and T-spaces is the same, 
so it is possible to couple them to N = I + T. Our problem then reduces to that of the case without strangeness: for- 
given / and S we have several possible values of N, for each of these we determine the linear combinations Eq. ( |54| ) 
with / replaced by N, and finally we uncouple / and T. We obtain 

\JJ Z ,II Z ,S}= ]T a J j? N3 {II 3 TT 3 \NN 3 }DL J3 Dl h Kl (58) 

J3N3I3T3 



3 Note that T = \S\/2. See below. 



where (II 3 TT 3 \NN 3 ) are the SU(2) Clebsch-Gordan coefficients. 

There is a further restriction of the possible collective states. Given a certain value of the baryon number B and 
the strangeness S, not all the values of isospin I are allowed. As discussed in Appendix B of Ref. [19], physical states 
should have hypercharge and isospin given by 

Y = B + S/3 = P -±^ ; I=\ (59) 

where p and q should be non- negative integer numbers. The allowed values of isospin / for states with S = 0, — 1 and 



— B are given in Table III, together with the corresponding values of T. Such values are obtained by imposing that 
the kaon wave function has to be completely symmetric under individual kaon exchange. 

It should also be noted that in the construction of the projector Eq. ( pq ) all the operations of G have to be taken 
into account (i.e. not only those of the proper subgroup). For this purpose the representations of the parity operation 
are also needed. For each baryon number, they are given in Table [II. Another important comment is that for odd 
baryon numbers the J and N quantum numbers are half-integers. For those cases one has to deal with the double 
group of G. 

V. NUMERICAL RESULTS 

In our numerical calculations we will use two standard sets of values for the Skyrme model parameters f^, e and m^. 
SET A corresponds to f n = 64.5 MeV, e = 5.45, m, = while SET Bto/ T = 54 MeV, e = 4.84, m^ = 138MeV @. 
In both cases we set the ratio fx/f-Tr to its empirical ratio Jk/Itt = 1-22. With these values we can calculate M so i, 
the kaon eigenenergies e„ and the radial integrals mi, m-i, d\ and di which appear in the expression of the moments 
of inertia and hyperfine splitting constants. The results are tabulated in Table IV. Using these values together with 
those for the angular integrals given in Table I all the parameters appearing in the collective Hamiltonians can be 
evaluated. For B = 1 we find that 9 = 1.01 fm and c = 0.50 for Set A and 9 = l.Olfm and c = 0.39 for Set B 
which provide a quite accurate description of the octet and decuplet baryon spectra plLB9[ . The numerical values 
of the parameters in the B = 2 collective Hamiltonian Eq. ( |47| ) are given in Table V. It interesting to compare the 
values of the inertia parameters with those obtained using the numerically obtained exact axially symmetric B = 2 
skyrmion H. For example, the corresponding values for Set B are 

K( = 30 M eV , K( = A8MeV, &i = 1.45 fm . (60) 

As we see the differences with the values listed in Table V are of only a few percent. On the other hand there not 
exist, so far, any calculation of hyperfine splitting constants using the exact numerical B = 2 skyrmion. Nevertheless, 
we can compare our results with those from a calculation based on an improved variational ansatz jl5[ which are, for 
Set B, 

c[ = 0.334 , c{ = 0.554 . (61) 

These values are also very similar to ours. This is also true for Set A. Taking into account that the corresponding 
inertia parameters are also very close to those given in Table V, it follows that our predicted dibaryon spectra coincide 
basically with the ones described in Ref. fla] . 

Results for the 5 = 3 — 9 inertia parameters and hyperfine splitting constants are listed in Tables VI and VII, 
respectively. As expected, the inertia parameters decrease with increasing baryon number. However, the decrease of 
the spin inertia appears to be much faster than that of the isospin one. This can be understood in the following way. 
Since we are interested in the overall behavior of inertias as a function of B we define, for both spin and isospin, the 
average value K = 1/3 ^ a K aa . As it can be seen from Table IV, while mi is roughly proportional to the baryon 
number, n%2 is basically independent of B. Therefore, assuming K w 1/9 and using Eqs.(jll| |l^) we have 

l/K J man TrC + b TrC , 1/K 1 man TrA + b TrA , (62) 

where a and b are constants roughly independent of n. On the other hand, it is not difficult to prove that the traces 
of the angular integrals appearing in these relations are given by 

TrA = 1 , TrA = 2n , TrC = 2n , TrC = 4X . (63) 

As shown in Ref. pOl, X < n 2 . In fact, from Table I we see that T is basically proportional to n 2 . Therefore, replacing 
Eq.(p3[) in Eq.(p2[) we obtain that K J should decrease as n 2 while K 1 goes only like 1/n. This behavior of the inertia 



parameters has important consequences in the multibaryon spectra. Namely, as the baryon number increases low 
lying non-strange states are expected to have the lowest possible value of isospin. For strange multibaryons this is 
not necessarily the case due to the coupling of the isospin to the kaonic spin T. 

The rotational energies for the non-strange multibaryons are given in Table VIII while those for S = — 1 states are 
given in Table IX and the corresponding to zero-hypercharge states in Table X. In all the cases we have included in 
the tables the lowest lying state and the first two excited states for each channel. Some general observations can be 
made. Due to the overall decrease of the inertia parameters the energy splittings become smaller as B increase. We 
also note that the ordering of the S = states is the same for both sets of parameters. For the S — — 1 states there 
is, however, one exception which corresponds to the second excited multibaryon with (B,S) = (6,-1). For Set A 
the second excited state is a 3 + while for Set B is a 2 + . It should be noticed, however, that the third excited states 
(not listed in Table IX) are precisely a 2 + for Set A and 3 + for Set B and that the energy difference with the second 
excited state is lMeV in both cases. For the Y — states the situation becomes more complicated as B increases. 
This is due to the rather small energy splittings between the different states. As a general trend we also note that 
the rotational energies are slightly smaller for Set B. This can be traced back to the fact that the moments of inertia 
are smaller for that set of parameters. 

As discussed above, for non-strange baryons the lowest lying states always have the lowest possible value of isospin. 
The corresponding spins are then given by the lowest value allowed by the symmetry constraints. As remarked in 
Ref. pQ| these values turn out to be consistent with those known for light nuclei with the exception of the odd 
values B = 5, 7,9. It should be stressed that at this point there is no obvious way to identify these rather compact 
multiskyrmion configurations with normal nuclei. Indeed, even for the B — 2 case it is not clear to which extent 
the deuteron wave function in the Skyrme model is represented by the torus configuration. Some analysis in terms 
of classical periodic orbits indicate that the two skyrmions spend most of their time at large separation and only a 
short time near the torus BG] . As the strangeness increases (in absolute value) the quantum numbers of the low- lying 
states become less obvious. This is a consequence of the interplay between the different terms in the corresponding 
collective Hamiltonian for non-zero values of T, In fact, the quantum numbers of the Y = states listed in Table X 
could be determined only after the calculation of the energies of a rather large set of allowed states. 

We discuss now the issue of the stability of the Y — I — states that we generically call multilambda states. The 



possible stability of a tetralambda state was first suggested in Ref. 1 31 1 . A similar conclusion was reached in Ref. |19f | 
where the existence of a stable heptalambda was also proposed. As already mentioned in the Introduction, in that 
work non-adiabatic corrections were neglected. We are now in position to check whether these effects do or do not 
affect the stability of these states. From Table X we observe that for Set B the g.s. Y = tetrabaryon is indeed a 
tetralambda state. This differs from the situation for Set A where the tetralambda is the first excited state. In any 
case, this does not affect the rotational contribution to the 4A— 2A mass difference. Using the energies given in Table 
X together with the values given in Table V for the parameters of H^ 2 (see Eq.(|47|)) we find that the rotational 
corrections decrease the binding by 36 MeV for Set A and by 26 MeV for Set B. These values are significantly smaller 
than the binding energy w 176 MeV obtained for both sets of parameters in the adiabatic approximation |]19). Thus, 
although the rotational corrections tend to decrease the binding, the tetralambda still turns out to be bound within 
the present approach. For the heptalambda we consider first its stability with respect to the decay into 3A + 4A. The 
non-adiabatic value of the corresponding binding energy is —177 MeV |l9(| . It should be noticed that the rotational 
energy of the zero-isospin (B, S) = (7, —7) state does not appear in Table X. In fact, the lowest lying of such states 
has E rot — 104MeV for Set A and E rot = 61MeV for Set B. That is, it shows up as an excited state with higher 
energy. Nevertheless, taking into account the rather large rotational energies of the Y = I = states with B = 3 
and 4 it happens that the binding energy of the heptalambda is increased by 45 MeV for both sets of parameters. 
For the case of the heptalambda ionization energy one can verify that the values given in Ref. |l9f remain basically 
unaffected by the rotational corrections. For this purpose one has to use the values of the rotational energies of the 
lowest Y = I = with B = 6. Such values (which are not listed in Table X) are 87 MeV for Set A and 52 MeV for 
Set B. 

VI. CONCLUSIONS 

In this work we have studied the non-adiabatic corrections to the masses of the multibaryons within the bound state 
approach to the SU(3) Skyrme model. To describe the multiskyrmion backgrounds we have used ansiitze based on 
rational maps. Such configurations are known to provide a good approximation to the exact numerical ones, and lead 
to a great simplification in the solution of the kaon eigenvalue equation. An important property of these approximate 
configurations is that they have the same symmetries as the exact ones. Consequently, the collective Hamiltonians 
and wave functions determined in this work are valid also in that case. They have been obtained making extensive use 
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of the properties of the corresponding symmetry groups. In particular, we have shown how the Finkelstein- Rubinstein 
phases fix, in a unique way, the one dimensional irreducible representations in which each wave function should 
transform. 

Using two standard sets of parameters for the effective SU(3) Skyrme action we have calculated all the inertia 
parameters and hypcrfine splitting constants for B < 9. We have found that as a general trend the isospin moments 
of inertia increase as n 2 while the spin ones as n, where B = n. Thus, the low lying non-strange multibaryons have 
the lowest possible value of isospin. The situation is more complicated in the case of strange particles for which there 
is a quite delicate interplay between the different terms contributing to the rotational energies. 

We have also estimated the non-adiabatic corrections to the tetralambda and heptalambda binding energies given 
in Ref. fl9[| . We found that these corrections are relatively small and do not affect the stability of these particles. This 
statement can be certainly extended to the recent studies on the stability of heavier flavored multiskyrmions [ p2| . 

We finish with a comment on the Casimir corrections to the multibaryon masses. Although these corrections are 
not expected to affect in any significant way the rotational energies obtained in the present work they might play 
some role in the determination of the multibaryon binding energies. Within the SU(2) Skyrme model it has been 
shown |33|] that they are responsible for the reduction of the otherwise large B = 1 soliton mass to a reasonable 
value when the empirical value of /^ is used. Here, we have avoided the B — 1 large mass problem by using the 
customary method of fitting f n to reproduce the nucleon mass |28[| . A more consistent approach should certainly use 
the empirical / w and include the Casimir corrections. In this respect, there have been recently some efforts [ p4[ to 
evaluate the corrections to the B = 1 mass in the SU(3) Skyrme model. Unfortunately, even in the SU(2) sector, 
almost nothing is known for B > 1. This is, of course, a very difficult task. Already in the SU(2) model, it requires 
the knowledge of the pion excitation spectrum around the non-trivial multiskyrmion up to rather large values of 
angular momentum. Nevertheless, recent studies of the SU(2) multiskyrmion low lying vibrational spectra p5| could 
be considered as first steps in this direction. 
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APPENDIX: 

In this Appendix we present the explicit expression of the rational maps used in this work. They are 

Ri = z , (Al) 

R 2 = z 2 , (A2) 

*3 = ^ ^ "J , (A3) 

z(z 2 - n/3) 

l + 2^^/3z 2 + z 4 

R 4 = = , (A4) 

1 - 2zV3 z 2 + z 4 V ' 

z (z ~ ib<sz 2 — arS) 

a^z^ + ib^z A — 1 

Re = z 2 ( t a 6 z 4 + 1) ' (A6) 

„ Z 5 — fl7 , . , 

R ^ = iwTT) ' (A7) 

R& = ~YT- 6 TT ' (A8) 

z z (iasz b — 1) 

z 3 (-z 6 + 3iV3z 4 + 9z 2 + 5iV3) + a 9 z (-iy/Zz 6 - z 4 + i^/iz 2 + l) 

Rc\ = = = ; = =- . (A9) 

5iV3z 6 + 9z 4 + 3iV3z 2 - 1 + a 9 z 2 (z 6 + iV3z 4 - z 2 - iy/3) 
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The numerical values of the real constants a%, bi appearing in these expressions are 

a 5 = 3.07, a 6 = 0.158, a 7 = 0.143 , a 8 = 0.137 , a 9 = 1.98 , b 5 = 3.94 . (A10) 

The reader can check that in most cases our maps agree with those given in Ref. [£0|. There are a few exceptions, 
however. For B = 7 we have choose a different orientation in the spin and isospin spaces in such a way that one of 
the 5-fold axes coincides with the z-direction. In the case of B — 9 we have selected the map for which the Td group 
operations are realized in exactly the same way in both spin and isospin spaces (namely, g — D g ). This is not the 
case for the B — 9 map given in Ref. [E0|. 
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TABLE I. Values of the diagonal elements of the angular integrals appearing in Eqs.(t 
listed in Table 1 are the values of 1 


$2j). Also 


n J 


A 


A 


B 


B 


C 


C 


1 1 


1/3 


2/3 


-2/3 


-4/3 


2/3 


4/3 


2 7T + 8/3 


I-tt/4 
I-tt/4 
it/2- 1 


4/3 





TV — 4 






-16/3 


7T-2 
7T-2 

8-2tt 


2tt 

2tt 

32/3 


3 13.58 


1/3 


2 


0.348 


4 


2 


18.11 


4 20.65 


0.391 
0.391 
0.219 


8/3 








8/3 


27.53 


5 35.75 


0.280 
0.280 
0.440 


10/3 


-0.090 
-0.090 
-0.051 


-3.649 
-3.649 
-2.038 


3.440 
3.440 
3.119 


50.65 
50.65 
41.71 


6 50.76 


0.356 
0.356 
0.285 


4 






0.089 






6.378 


4.137 
4.137 
3.725 


71.82 
71.82 
59.40 


7 60.87 


1/3 


14/3 








14/3 


81.16 


8 85.63 


0.312 
0.312 
0.376 


16/3 






-0.044 






-8.91 


5.171 
5.171 
5.658 


105.89 
105.89 
130.74 


9 113.07 


1/3 


6 


-0.062 


-7.51 


6 


150.76 



TABLE II. Symmetry group G, generators of the proper subgroup, their corresponding FR 
phases and the parity operations for B=3-9. The directions of the 3-fold axes in B — 7 

are defined by the spherical angles (4> a ,0 a ) = ( tt/5, 
I 37r/5, arccos 



V(5 + 2V5)/15 



and (<j>p,dp) 



1/^15 + 6^5 



B 



G 



Generators of proper subgroup and FR phases 



{9i,D 91 } 



lax 



{92,D g2 } 



792 



Parity 
operation 



3 


T d 


(fixyz y-ixyz-i 


1 


{CI, CI} 


4 


o h 


{(%»,(%} 


1 


{CI, CI} 


5 


D 2d 


{CI, CI} 


1 


{CI, CI} 


6 


D 4d 


{C%,C%} 


-1 


{cv,ci} 


7 


h 


{CI, (CI) 3 } 


1 


{CS,(C§) a } 


8 


D 6d 


{CI, CI} 


1 


{CI, CI} 


9 


T d 


(fjxyz rjxyz-. 


1 


{CI, CI} 



1 


{CI, CI} 


1 


{E,CI} 


-1 


{C 4 , C 4 } 


-1 


{CI, CI} 


1 


{E,E} 


1 


{Cl2, Cg} 


1 


{CI, CI} 
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TABLE III. Allowed values of I and T for states with dif- 
ferent strangeness for B=3-9. 



B 


S 


I 


T 


3 





1/2, 3/2,..., 9/2 







-1 


0, 1,...,4 


1/2 




-3 


0, 1,...,3 


3/2 


4 





0,1,..., 6 







-1 


1/2, 3/2,..., 11/2 


1/2 




-4 


0, 1,...,4 


2 


5 





1/2, 3/2,..., 15/2 







-1 


0, 1,...,7 


1/2 




-5 


0, 1,...,5 


5/2 


6 





0,1,..., 9 







-1 


1/2, 3/2,..., 17/2 


1/2 




-6 


0, 1,...,6 


3 


7 





1/2, 3/2,..., 21/2 







-1 


0, 1,. ..,10 


1/2 




-7 


0,1,..., 7 


7/2 


8 





0, 1,...,12 







-1 


1/2, 3/2,. . . ,23/2 


1/2 




-8 


0, 1,...,8 


4 


9 





1/2, 3/2,. . . ,27/2 







-1 


0, 1,. ..,13 


1/2 




-9 


0, 1,...,9 


9/2 



TABLE IV. Numerical values of the radial integrals appearing in the expressions of moments of 
inertia and hyperfine splitting constants. 







SETA 






SETB 




B 


mi(/m) 


m 2 (fm) 


rfi 


d 2 (x .01) 


mi(fm) 


m 2 (fm) 


di 


d 2 (x .01) 


1 


1.27 


0.212 


0.227 


2.29 


1.22 


0.304 


0.275 


2.93 


2 


1.95 


0.233 


0.233 


1.27 


2.21 


0.335 


0.272 


1.61 


3 


2.58 


0.241 


0.228 


0.89 


3.06 


0.349 


0.261 


1.13 


4 


3.00 


0.246 


0.213 


0.69 


3.65 


0.359 


0.243 


0.89 


5 


3.74 


0.249 


0.217 


0.57 


4.53 


0.363 


0.244 


0.73 


6 


4.32 


0.251 


0.214 


0.48 


5.23 


0.366 


0.239 


0.62 


7 


4.65 


0.253 


0.205 


0.43 


5.67 


0.371 


0.229 


0.55 


8 


5.39 


0.254 


0.208 


0.38 


6.51 


0.371 


0.230 


0.49 


9 


6.09 


0.255 


0.209 


0.34 


7.30 


0.371 


0.230 


0.44 



TABLE V. Parameters for B=2 



SET 


K( {MeV) 


K{ (MeV) 


e| (fm) 


c[ 


d 


A 
B 


33.42 
27.63 


53.68 
45.20 


1.15 
1.40 


0.409 
0.306 


0.631 
0.562 
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TABLE VI. Inertia parameters for B=3-9 







SET A 






SET B 




B 


K J {MeV) 


K 1 {MeV) 


K m {MeV) 


K J {MeV) 


K 1 {MeV) 


K m {MeV) 


3 


15.23 


50.77 


9.55 


12.11 


41.03 


7.80 


4 


8.66 


39.70 
39.70 
32.88 





6.72 


30.98 
30.98 
25.89 





5 


5.20 


28.29 


-1.17 


4.03 


22.30 


-0.96 




5.20 


28.29 


-1.17 


4.03 


22.30 


-0.96 




5.88 


33.91 


-0.89 


4.57 


26.47 


-0.73 


6 


3.67 


26.12 





2.84 


20.45 







3.67 


26.12 





2.84 


20.45 







4.25 


24.48 


1.23 


3.31 


19.28 


1.04 


7 


3.09 


23.06 





2.38 


17.90 





8 


2.39 


19.48 





1.85 


15.28 







2.39 


19.48 





1.85 


15.28 







2.11 


21.08 


-0.61 


1.63 


16.50 


-0.52 


9 


1.78 


17.75 


-0.43 


1.39 


14.02 


-0.36 



TABLE VII. Hyperfine splitting constants for B=3-9 





SETA 




SETB 




B 


c J 


c' 


c J 


c> 


3 


-0.62 


0.55 


-0.64 


0.48 


4 





0.55 
0.55 
0.46 





0.48 
0.48 
0.37 


5 


0.22 


0.48 


0.23 


0.41 




0.22 


0.48 


0.23 


0.41 




0.15 


0.57 


0.16 


0.51 


6 





0.53 





0.46 







0.53 





0.46 




-0.28 


0.49 


-0.30 


0.43 


7 





0.53 





0.46 


8 





0.51 





0.45 







0.51 





0.45 




0.28 


0.55 


0.31 


0.49 


9 


0.23 


0.52 


0.25 


0.46 



15 



TABLE VIII. Quantum numbers and rotational energies 
for S — states 







SET A 






SETB 




B 


J F 


I 


N E rot [MeV] 


J F 


/ 


N E rot [MeV] 


3 


1/2+ 


1/2 


1/2 


64 


1/2+ 


1/2 


1/2 


52 




5/2~ 


1/2 


1/2 


147 


5/2~ 


1/2 


1/2 


117 




3/2" 


3/2 


3/2 


205 


3/2" 


3/2 


3/2 


164 


4 


0+ 











0+ 













4+ 








173 


4+ 








134 




0+ 


2 


2 


238 


0+ 


2 


2 


186 


5 


1/2+ 


1/2 


1/2 


28 


1/2+ 


1/2 


1/2 


22 




3/2+ 


1/2 


1/2 


40 


3/2+ 


1/2 


1/2 


31 




3/2" 


1/2 


1/2 


44 


3/2~ 


1/2 


1/2 


34 


6 


1+ 








7 


1+ 








6 




3+ 








41 


3+ 








34 




0+ 


1 


1 


52 


0+ 


1 


1 


41 


7 


7/2+ 


1/2 


1/2 


66 


7/2 + 


1/2 


1/2 


51 




3/2+ 


3/2 


3/2 


98 


3/2+ 


3/2 


3/2 


76 




9/2+ 


3/2 


3/2 


163 


9/2+ 


3/2 


3/2 


126 


8 


0+ 











0+ 













2+ 








14 


2+ 








11 




1+ 


1 


1 


44 


1+ 


1 


1 


34 


9 


1/2+ 


1/2 


1/2 


14 


1/2+ 


1/2 


1/2 


11 




5/2- 


1/2 


1/2 


30 


5/2- 


1/2 


1/2 


24 




7/2- 


1/2 


1/2 


39 


7/2- 


1/2 


1/2 


31 



TABLE IX. Quantum numbers and rotational energies for S — — 1 states 









SETA 








SET B 




B 


J" 


I 


N 


E rot [MeV] 


r 


/ 


N 


E rot [MeV] 


3 


1/2+ 





1/2 


38 


l/2 + 





1/2 


29 




1/2+ 


1 


1/2 


84 


1/2+ 


1 


1/2 


72 




5/2- 





1/2 


122 


5/2" 





3/2 


95 


4 


0+ 


1/2 


1/2 


6 


+ 


1/2 


1/2 


7 




4+ 


1/2 


1/2 


180 


4+ 


1/2 


1/2 


141 




0+ 


3/2 


3/2 


191 


0+ 


3/2 


3/2 


144 


5 


1/2+ 





1/2 


11 


1/2+ 





1/2 


7 




3/2+ 





1/2 


23 


3/2+ 





1/2 


17 




3/2- 





1/2 


28 


3/2- 





1/2 


21 


6 


1+ 


1/2 





12 


1 + 


1/2 





10 




0+ 


1/2 


1 


32 


0+ 


1/2 


1 


23 




3+ 


1/2 





49 


2+ 


1/2 


1 


37 


7 


7/2+ 





1/2 


54 


7/2 + 





1/2 


40 




3/2+ 


1 


1/2 


74 


3/2+ 


1 


1/2 


56 




7/2+ 


1 


1/2 


75 


7/2+ 


1 


1/2 


60 


8 


0+ 


1/2 





3 


+ 


1/2 





3 




2+ 


1/2 





18 


2+ 


1/2 





14 




1+ 


1/2 


1 


28 


1+ 


1/2 


1 


21 


9 


1/2+ 





1/2 


4 


1/2+ 





1/2 


3 




5/2- 





1/2 


20 


5/2- 





1/2 


15 




1/2+ 


1 


1/2 


21 


1/2+ 


1 


1/2 


18 



1G 



TABLE X. Quantum numbers and rotational energies for Y = states 









SETA 








SET B 




B 


J r 


I 


N 


£ ro *[MeV] 


r 


/ 


N 


£ rot [Mel/] 


3 


1/2+ 


1 


1/2 


50 


1/2+ 


1 


1/2 


45 




3/2~ 





3/2 


77 


3/2~ 





3/2 


52 




3/2 - 


1 


3/2 


123 


5/2+ 





3/2 


89 


4 


+ 


2 


2 


51 


0+ 





2 


43 




0+ 








72 


0+ 


2 





54 




0+ 


1 


1 


109 


0+ 


1 


2 


77 


5 


l/2 + 


1 


3/2 


29 


l/2 + 


1 


3/2 


21 




1/2" 


1 


3/2 


32 


1/2" 


1 


3/2 


23 




1/2+ 


2 


1/2 


39 


3/2+ 


1 


3/2 


30 


6 


+ 


2 


1 


24 


0" 


1 


2 


16 




o- 


1 


2 


26 


1- 


1 


2 


22 




1- 


1 


2 


33 


1+ 


1 


2 


23 


7 


3/2+ 


2 


7/2 


32 


3/2+ 


2 


7/2 


28 




5/2+ 


1 


7/2 


65 


5/2+ 


1 


7/2 


42 




7/2+ 


1 


7/2 


87 


7/2+ 


1 


7/2 


59 


8 


+ 


2 


2 


19 


0+ 


2 


2 


16 




2+ 


2 


2 


31 


2+ 


2 


2 


24 




2+ 


2 


2 


33 


2+ 


2 


2 


27 


9 


1/2" 


2 


5/2 


25 


1/2" 


2 


5/2 


18 




3/2" 


2 


5/2 


29 


3/2" 


2 


5/2 


21 




3/2+ 


3 


3/2 


31 


3/2+ 


2 


5/2 


24 



17 



